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ABSTRACT 

PKS 0558-504 is the brightest radio-loud Narrow-Line Seyfert 1 galaxy at X-ray energies. Here we 
present results from the radio, optical, UV, and X-ray bands obtained with SWIFT, XMM-Newton, 
and ATCA during a 10-day monitoring campaign in September 2008. The simultaneous coverage 
at several wavelengths makes it possible to investigate in detail the broadband spectral energy dis- 
tribution (SED) and the energetic of this source. The main results can be summarized as follows. 
The ATCA reveals the presence of an extended radio emission in PKS 0558-504 with two lobe- like 
structures at ^7" from the bright central source. The extended radio structure and the low value of 
the radio-loudness similar to radio-quiet Seyfert galaxies coupled with constraints from higher energy 
bands argue against a jet-dominated emission in PKS 0558-504. The study of the SED, which is 
dominated by a nearly constant optical-UV emission, supports the conclusion that PKS 0558-504 is 
accreting at super-Eddington rate. This conclusion was reached assuming Mbh = 2.5 x 10* Mq, which 
was obtained with a new scaling method based on X-ray spectral variability results. A comparison 
between the accretion luminosity and the kinetic power associated with the jet suggests that in this 
source the accretion power dominates in agreement with the results obtained from Radiation-MHD 
simulations of Galactic black holes (GBHs) accreting at the Eddington rate. The combined find- 
ings from this panchromatic investigation strongly suggest that PKS 0558-504 is a large scale analog 
of GBHs in their highly accreting intermediate state. Importantly, PKS 0558-504 may also be the 
prototype of the parent population of the very radio-loud NLSls recently detected at 7-ray energies. 

Subject headings: Galaxies: active - Galaxies: jets - Galaxies: nuclei - X-rays: galaxies 



1. INTRODUCTION 

Bipolar relativistic jets are common features in a va- 
riety of astrophysical objects, most notably in Galactic 
Black Holes (GBHs) and Active Galactic Nuclei (AGNs). 
Accretion of gas onto black holes is thought to power 
these coUimated outflows. However, the details of the 
jet formation as well as the nature of the coupling be- 
tween accreting matter and outflows are still among the 
outstanding open questions in high energy astrophysics. 

Because of their vicinity and hence their high bright- 
ness, the temporal and spectral properties of GBHs are 
much better known and can be used to infer informa- 
tion on their more powerful, extragalactic analogs. In- 
deed considerable progress in this field has been made 
by multi-wavelength correlated studies of GBHs in dif- 
ferent spectral states, which allow one to study the 
link between accretion (generally probed by the X-ray 



emission) and jet properties (in the radio regime) on 
"human" timescales. It is now well established that 
GBHs undergo state transitions, switching between two 
main states: the low/hard (LS) and the high/soft state 
(HS) passing through soft and hard intermediate states 
(IS), which are also named very high states (VHS) 
when they occur at high values of accret i on ra te (see 
iMcCIintock fc RemiUardl l2006t iDone et all l2007l for re- 
cent comprehensive reviews on GBHs). Each spectral 
state is unambiguously characterized by a specific com- 
bination of X-ray temporal and s pectral properties a nd 
by well defined radio features (see lFender et al.ll2004D . 

In the study of the disk-jet link one of the most inter- 
esting spectral states is the IS/VHS, which is generally 
characterized by powerful transient relativistic ejections 
in the radio coupled with highly variable X-ray emission 
that is unambiguously associated with the accretion fiow, 
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unlike the LS where the origin of the X-rays is still matte r 
of debate (e.g. lMarkoff et al.|[200l IZdziarski et al.ll2004[ ). 

Unfortunately, the physical conditions that lead to the 
transient relativistic ejections during the IS/VHS are still 
poorly understood mostly because of their short duration 
in GBHs. Since the dynamical time scales are propor- 
tional to the black hole mass, in individual AGNs it is 
not possible to observe long-term phenomena occurring 
in GBHs such as the canonical spectral transitions. On 
the other hand, AGNs may provide better constraints on 
short-lived GBH phenomena and hence shed light on the 
jet formation and the interplay between accretion and 
ejection processes. 

In the framework of the AGN-GBH unification, the 
Narrow-Line Seyfert 1 galaxies (NLSls) are the best 
candidates for large-scale analogs of GBHs in the 
IS/VHS. NLSls are historically identified by their op- 
tical emission line properties: the ratio [O IIIJ/H/3 is 
less than 3 and FW HM H/3 is less than 2000 km s'^ 
(jOsterbrock fc^Q ggc 1985; Goodrich 1989). They are 
seldom radio loud (JKomossa et al.ll2006[ ). although recent 
studies reveal the existence of severa l NLSls character - 
ized by very high radio-loudness (e.g. lYuan et all 120081 ). 
Recently, a few of these very radio-loud NLSls have been 
detected at 7-ray energies by the feryni/L^T collabora- 
tion confirming that these sources posses relativistic jets 
observed at sinall vi ewing angles (jAbdo et al.ll2009al lbll3: 
iFoschini et al.ll2009D . In the X-rays, NLSl are gener- 
ally characterized by steep spectra and strong va riabil- 
ity (e .g. Brandt et al. 1999; Lcighly 1999a b; Gruo e et al.l 
[200 It ). Based on these properties, it has be en suggested 
that NLSl are AGN in their early phase (jGrupe et al.l 
II999I : iMathud l2000( ) , characterized by relati vely small 
black hole masses (e.g. lGrupe fc Mathur" 200 J ) . and very 
high accretion rates in terms o f Eddingto n units (e.g. 
Borqson fc Greenlll992l:ISule ntic ct al. 2000; Grup e_et al.l 



2010t see also lMarconi et al.ii2008: .Decarh et al...2008l for 



a discording view). 

PKS 0558-504 {z = 0.137) is the brightes t radio-loud 
{R = i^(6 cm)/L^(B) =i 27. iSiebert etall [l999) NLSl 
in the X-ray band and therefore one of the best stud- 
ied at these energies. Previous X-ray observations with 
different satellites have confirmed that PKS 0558-504 
shows the characteristic NLSl properties: strong vari- 
ability, steep X-ray spectrum, substantial soft excess, and 
relatively high luminosity. By comparing the X-ray ob- 
servations from several satellites over a decade, it is evi- 
dent that on long timescales the strong X-ray variability 
of PKS 0558-504 occurs persistently but without being 
accompanied by significant spectral variability: the spec- 
trum above 2 keV is consistently described by a power- 
law m odel with photon index F ~ 2.2 (G liozzi et al.l 
l2000t) . This conclusion has been confirmed by a recent 
monitoring campaign with RXTE: the long-term achro- 
matic variability appears to be consistent with the behav- 
ior of GBHs during the transition LS-to-HS and inconsis- 
tent wi th the spectral varia bility typical of jet-dominated 
AGNs (|Gliozziet al.l[2007l) . 

In the last ten years, PKS 055 8-504 has been ob- 
served repeatedly by XMM-Newton dO'Brien et al.|[200ll: 
iGliozzi et ahl 120011; iBrinkmann et al.l I2004D . The most 
recent and highest quality data, obtained with the 
XMM-Newton EPIC cameras from five consecutive or- 
bits in September 2008, confirm that PKS 0558-504 is 



highly variable on all sampled timescales, and a spec- 
tral analysis indicates that the 2-10 keV energy range 
is well fitted with a power law whose slope steepens 
as the source brightens and the soft excess is well de- 
scribed by a low-temp erature Comptonization model 
(jPapadakis et al.ir2010al ). These findings suggest that de- 
spite its radio loudness PKS 0558-504 appears to behave 
like a "normal" Seyfert galaxy from the X-ray point of 
view. 

The main goal of this work is to shed light on the 
energetic of this powerful source using multi-wavelength 
data. For this purpose, after describing the observations 
and data reduction in § 2, we first try to assess the role 
played by the jet emission beyond the radio band (in § 3) 
and then to constrain the black hole mass in PKS 0558- 
504 without the use of optical measurements (in § 4). In 
Section 5 we describe the multi-wavelength behavior of 
PKS 0558-504 based on a 10-day multi- wavelength cam- 
paign carried out with the SWIFT XRT and UVOT si- 
multaneously with the deep XMM-Newton observations 
in September 2008, and complemented with one radio 
observation with the Australia Telescope Compact Ar- 
ray (ATCA). Finally, in § 6 we discuss the implications 
of the main results and draw our conclusions. Hereafter, 
we adopt a cosmology with Hq = 71 km s~^ M pc~^, 
r^A = 0.73 and ^m = 0.27 ()Bennet et al.|[200l ; with 
the assumed cosmological parameters, the luminosity dis- 
tance of PKS 0558-504 is 642 Mpc. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. SWIFT 

PKS 0558-504 was obser ved by the SWIFT Ga.mma.- 
Ray Burst explorer mission (jGehrels et al.|[200^ between 
2008 September 9 - 16. The details of this intense mon- 
itoring campaign are summarized in Tables 1 and 2 . 
The SWIFTX-ia.y Telescope fXRT: lBurrows et al.l[2005l) 
observations were all perfor med in Windowed Timing 
mode (WT; iHiU et al.l [200l in order to avoid the ef- 
fects of pile-up. Data were reduced by the task xrt- 
pipeline version 0.12.1., which is included in the HEA- 
SOFT package 6.8. Source and background photons were 
selected in boxes 40 pixel long. Only single and double 
events (GRADES to 2) were used. Source photons 
for the light curve and spectra were extracted with XS- 
ELECT. The auxiliary response files (ARFs) were cre- 
ated using xrtmkarf version 0.5.6 and the response matrix 
swxwt0to2s6_2001 01 01 vOl 1 . rmf. 

The UV/Optical Telescope (UVOT; Roming et al. 
2005) observed PKS 0558-504 in aU 6 fihers. The ex- 
posure times in each of the filters per segment are given 
in Table 1 and the corresponding magnitudes in Table 
2. The U VOT data were re duced and analyzed as de- 
scribed in iPoole et al.l (|2008D . Source photons were ex- 
tracted from the coadded images files for each segment 
with a radius of 5" Backgound photons were selected in 
a nearby source-free region with r=20'.' Magnitudes and 
fiuxes were measured using the UVOT tool uvotsource. 
All magnitudes listed in Table 2 were corrected for Galac- 
tic reddening f^R_v=0.044; iSchlegel eTall [l998l ) . The 
correction factors we re based on the stand ard reddening 
correction curve s by iCardelli et al.l ( 1989[ ) as described 
by equation 2 in lRoming et al.l (|2009[ ). 



2.2. XMM-Newton 
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TABLE 1 
SWIFT Observation log of PKS 0558-504 





Segment 


Start time 


End Time 


M,ID 




Observing 


time given in s 






(UT) 


(UT) 




TxRT 


Tv 


Tb 


Tu 


Tuvwi 


TuVM2 


TuVW2 


001 


2008-09-07 08:10 


2008-09-07 11:36 


54716.41 


2215 


122 


122 


122 


244 


210 


488 


002 


2008-09-08 01:50 


2008-09-08 06:50 


54717.20 


2429 


199 


199 


199 


399 


358 


798 


003 


2008-09-09 09:52 


2008-09-09 11:48 


54718.45 


2327 


191 


191 


191 


381 


532 


763 


004 


2008-09-10 00:15 


2008-09-10 18:20 


54719.39 


2203 


181 


181 


181 


361 


509 


724 


005 


2008-09-11 13:21 


2008-09-11 19:59 


54720.69 


1897 


163 


163 


163 


326 


344 


653 


006 


2008-09-12 00:36 


2008-09-12 05:37 


54721.15 


2008 


155 


155 


155 


310 


483 


621 


007 


2008-09-13 15:09 


2008-09-13 18:40 


54722.68 


1869 


148 


148 


148 


296 


447 


594 


008 


2008-09-14 00:42 


2008-09-14 05:40 


54723.13 


2363 


240 


240 


240 


483 


682 


966 


009 


2008-09-15 16:56 


2008-09-15 21:54 


54724.80 


2353 


184 


184 


184 


369 


557 


739 


010 


2008-09-16 07:20 


2008-09-16 12:11 


54725.90 


2072 


165 


165 


165 


330 


461 


660 





PKS 0558-504 was observed by XMM-Newton from 
September 9, 2008, 1:19 UT to September 16, 2008, 12:02 
UT. The EPIC pn and MOSl cameras were operated 
in small window mode, the M0S2 in timing mode, and 
both RGS in spectroscopy mode. In the present work 
we only use EPIC pn data, which were processed with 
the XMM-Newton Science Analysis Software (SAS) 8.0. 
The recorded single and double events were screened to 
remove known hot pixels and other data flagged as bad: 
only data with PATTERN<4, FLAG=0 were used. A more 
detailed descrip tion of the XMM-Newto n data analysis 
can be found in iPapadakis et aLl (|2010al ). 

2.3. ATCA 

Radio monitoring observations of the PKS 0558-504 
with the ATCA commenced in April 2006. The observa- 
tion session that overlapped with the SWIFT a.nd XMM- 
Newton monitoring campaign was carried out on Septem- 
ber 16, 2008, 22:50 UT. The source showed a slightly 
inverted spectrum between 18.5 and 4.8 GHz (a ~ 0.3, 
where the flux density /^ is related to the spectral index 
a by the relation f^, oc v^°'). The peak flux density of 
0.09±0.01 Jy was measured at 4.8 GHz. Our previous 
images at this frequency indicate that almost 80% of ra- 
dio intensity originates from the bright unresolved core, 
as shown in Figure [TJ No flux density fluctuations were 
detected on intra-hourly timescales. 

2.4. SED fitting 

The broadband spectral analysis w as performed using 
the XSPEC V. 12.4 software package (|Arnaudlll996l ). We 
used FLX2XSP in FTQQLS to transform the optical and 
UV fluxes into suitable units for XSPEC. All the X-ray 
spectra were re-binned so that each bin contained at least 
20 counts for the x^ statistic to be valid. 

3. RADIO PROPERTIES OF PKS 0558-504 

Despite the fact that PKS 0558-504 has been one of 
the first radio-loud NLSl ever discovered, not until this 
project its radio properties and activity have been stud- 
ied in deserving detail. Our multi-wavelength, radio 
monitoring and imaging program with the ATCA and 
the VLBI shows the long-term variability pattern sim- 
ilar to t he flat spectrum, low luminosity compact radio 
sources (jFalcke et al.ll2000l) . A detailed description of the 
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Fig. 1. — The ATCA image with the rms of 2 X lO""* Jy/bcam was 
obtained on the 19th of January 2008. The superimposed contour lines 
show isophotes between 2% and 95% of the peak flux density. The 
beam size is shown in the bottom left corner. 



radio properties will be presented elsewhere (Kedziora- 
Chudczer et al. 2010, in preparation). Here we limit our- 
selves to briefly discuss the morphology of the radio emis- 
sion and the intensity of the source during the SWIFT 
and XMM-Newton campaigns in September 2008. 

From Figure [T] it is evident that the radio emission 
of PKS 0558-504 resembles an FRI radio galaxy that is 
characterized by two diffuse lobes located at ^ 7" from 
the bright, central core. At the redshift of z = 0.137 
the projected linear size of the full structure is close to 
46 kpc. Interestingly, the extended radio emission ob- 
served in PKS 0558-504 is at odds with the compact 
radio structure generally observed in other radio-loud 
NLSls, which are characte rized by radio-lou dness param- 
eters larger than 100 (e.g. lYuanet al.ll2008D . 

Historically, the radio-loudness parameter has been de- 
fined as i? = Li/(6 cm)/Lj^(B) with i? = 10 considered as 
the conventional bo undary between radio-lo ud and radio- 
quiet AGN s (e.g . iKellermann et all |1994[ ). However, 
iHo fc PengI (|2001t ) demonstrated that, when the nuclear 
luminosities are properly measured, most of the Seyfert 
galaxies appear to be radio-loud, indicating that i? = 10 
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TABLE 2 

SWIFT XRT COUNT RATES AND HARDNESS RATIOS AND UVOT MAGNITUDES OF PKS 0558-504 



Segment XRT rate XRT HR 



V 



B 



U 



UVWl 



UVM2 



UVW2 



001 
002 
003 
004 
005 
006 
007 
008 
009 
010 



0.747±0.020 
1.406±0.025 
0.933±0.021 
0.856±0.025 
0.808±0.022 
1.080±0.025 
1.010±0.025 
0.777±0.019 
0.577±0.017 
0.840±0.022 



-0.11±0.02 
-0.05±0.02 
-0.12±0.02 
-0.05±0.03 
-0.15±0.03 
-0.11±0.02 
-0.07±0.02 
-0.07±0.02 
-0.13±0.03 
-0.16±0.03 



14.90±0.03 
14.85±0.02 
14.98±0.02 
14.90±0.03 
14.84±0.02 
14.84±0.02 
14.88±0.02 
14.89±0.02 
14.86±0.02 
14.87±0.02 



15.06±0.02 
15.07±0.01 
15.08±0.01 
15.07±0.02 
15.05±0.01 
15.05±0.01 
15.02±0.01 
15.02±0.01 
15.05±0.01 
15.04±0.01 



13.83±0.01 
13.85±0.01 
13.85±0.01 
13.89±0.01 
13.84±0.01 
13.86±0.01 
13.85±0.01 
13.84±0.01 
13.83±0.01 
13.87±0.01 



13.54±0.01 
13.55±0.01 
13.63±0.01 
13.62±0.01 
13.57±0.01 
13.53±0.01 
13.54±0.01 
13.55±0.01 
13.56±0.01 
13.58±0.01 



13.23±0.02 
13.22±0.01 
13.36±0.01 
13.33±0.01 
13.23±0.01 
13.25±0.01 
13.24±0.01 
13.23±0.01 
13.20±0.01 
13.23±0.01 



13.29±0.01 
13.30±0.01 
13.43±0.01 
13.33±0.01 
13.33±0.01 
13.30±0.02 
13.30±0.01 
13.31±0.01 
13.30±0.01 
13.31±0.01 



The magnitude was corrected for reddening with i?B-v— 0.044 given bv lSchlegel et all ] 19981) . The errors given in this table are statistical errors. 



is not a reliable boundary. This conclusion is furthe r 
supported by the recent findings of iSikora et al.l ()2007[ ). 
who showed that radio-loud and radio-quiet AGNs form 
two distinct branches when the radio-loudness is plotted 
versus the Eddington ratio AEdd- Importantly, while the 
two branches are nearly horizontal and well separated 
at low values of the Eddington ratio, as AEdd increases 
the branch slopes become negative and the branches 
broaden with substantial overlap at large values of the 
Eddington ratio. Unfortunately, the super-Eddington 
accretion rate in PKS 0558-504 (see Sect. 5.2) cou- 
pled with the lack of a clear distinction between the two 
branches at large values of the Eddington ratio hampers 
the study of the radio-loudness for PKS 0558-504 using 
the logi? — log AEdd plot. 

Nevertheless, the radio-loudness of PKS 0558-504 can 
be qua ntitatively asse s sed by comparing it with the find- 
ings of iPanessa et al.l ()2007D who carried out a detailed 
investigation of the radio-loudness in two large samples 
of radio-quiet Seyfert galaxies and Low-Luminosity Ra- 
dio Galaxies (LLRGs). In this study both the classical 
radio loudness parameter R and the X-ray radio loudness 
i?x = J^ii/(6 cm)/L2-io keV wer e used; the latter par am- 
eter was introduced by Terash ima fc Wilsonl (|2003[ ) to 
circumvent extinction problems that usually affect the 
optical emission and may lead to overestimates of R. 
Using nearly simultaneous radio, optical, and X-ray ob- 
servations (we have used the flux values measured during 
the last day of the XMM-Newton - 5" FF/FT campaign be- 
cause they are almost contemporaneous with the ATCA 
observation) for the logR - logi?x plot, PKS 0558-504 
appears to be fully consistent with radio-quiet Seyfert 
galaxies and inconsistent with radio-loud objects, as vi- 
sually demonstrated in FiglH wher e the b oundaries have 
been determined bv IPanessa et al.l ()2007l ). 

In conclusion, the extended radio emission and the 
symmetric location of the lobes in PKS 0558-504, which 
are reminiscent of the typical structure observed in FRI 
radio galaxies, argue against a highly beamed source. 
This is confirmed by the simultaneous radio-loudness 
parameters that are in full agreement with radio-quiet 
Seyfert galaxies. On the other hand, our VLBI imaging 
at 2.3 GHz (Kedziora-Chudczer et al. 2010, in prepara- 
tion), which shows a one-sided pc-scale jet, suggests that 
the presence of beamed emission cannot be ruled out for 
PKS 0558-504 in the radio band. For the sake of clarity 
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Fig. 2. — Classical radio loudness R plotted versus the X-ray ra- 
dio loudness i?x. PKS 0558—504 is well within the radio-quiet region 
and apparently inconsistent with the radio-loud objects. The bound- 
aries between rad i o-loud and radio-quiet objects have been determined 
bv .Panessa et al.l 120071) using two large samples of radio-quiet Seyfert 
galaxies and LLRGs. 



and for historical reasons in the reminder of the paper 
we will maintain the classification of radio-loud NLSl for 
PKS 0558-504. 

4. BLACK HOLE MASS OF PKS 0558-504 

In order to determine the accretion rate of PKS 0558- 
504 and investigate its energetics, it is first necessary to 
constrain the mass of the supermassive black hole. Here, 
we report four different and independent measurements 
of A/bh for PKS 0558-504. 

• (1) As for many AGNs for which there are no direct 
measurements from reverberation mapping, the Mbh in 
PKS 0558-504 has been estimated from the virial re- 
lationship Mbh — f Ri^VY / G (where / is an unknown 
geometric factor, and R and (AV)^ are the radius and ve- 
locity dispersion of broad line r egion, respectively). Thi s 
yielded A/bh ~ 6 x lO'^ Afg (jPapadakis et al.ll2010aD . 
which is relatively small when c ompared to the typica l 
masses of radio- loud AGNs (e.g. iMclure fc Jarvisl[200^ . 
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but fully consistent with t he values derived in N LSls 
with the same method (e.g. iGrupe fc Mathuill2004f) . 
Given that the optical measurements -L 

(jCorbin fc SmithI [2000) and Hp (jCorbin I [l99l - are 
non-simultaneous, and given the uncertainty on the 
geometric factor / and the controversy about the 
applicatio n of the virial method to NLSls (s ee e.g. 
Marc oni e t al.l 120081 : iDecarh et all 120081 : but also lNetzed 
20091) . it is important to constrain A/bh also with 
alternative techniques that are independent of optical 
measurements and any assumption on the BLR. 

• (2) One possible alternative method is based on the 
so called "fundamental p l ane" of bla ck holes (BHs) intro - 
duced bv iMerloni et all (|2003D and iFalcke et al.l (|2004[ ). 
where Mbh is related to both the X-ray and radio lu- 
minosities in any accreting BH s ystems. Recent l y this 
relationship has been refined by iGiiltekin et al.l ( 2009f ) 
by utilizing only sources with BH masses that have been 
determined dynamically. Using the latter relationship 
in combination of quasi simultaneous measurements of 
the radio and X-ray emission in PKS 0558-504 we derive 
Mbh — 3 X 10* M©, which is larger than the value ob- 
tained from the virial theorem by a factor of ~4. This 
finding does not depend on the nearly simultaneous na- 
ture of the observations, since we derive a very similar 
value for Mbh by using radio and RXTE luminosities 
averaged over an interval of one year. 

• (3) For AGNs with evenly-sampled long-term cov- 
erage in the X-ray band, a viable technique to estimate 
the mass of the black hole is based on the relationship 
logMBH = (logTB-F0.981ogLboi + 2.32)/2.1, where Mbh 
is in 10^ Mq units, Lboi is the bolometric luminosity in 
units of 10*'* erg s~*, and Tb the time break in days de- 
rived from power spectral density (McHardy et al. 2006). 
PKS 0558-504 has been regularly monitored with RXTE 
since March 2004, making it possible to estimate Mbh 
with the above formula. The main findings of a detailed 
temporal study obtained combining the long-term RXTE 
light curve and the deep XMM-Ne wton obseryation in 
September 2008 are reported bv Papadakis et all ()2010b[ ) 
and suggest that Mbh ~ 2 - 3 x 10* Mq. 

• (4) An additional independent way to determine 
the mass in BH systems relies on the fact that hard 
X-rays are produced by the Comptonization process in 
both stellar and supermassive b lack holes. Recently, 
IShaposhnikov fc Titarcbli] (|2009t ) discovered that GBHs 
present a universal scalable relationship between the 
photon index and the normalization of the Bulk Mo- 
tion Comptonization (BMC) model during their spectral 
transitions. They also demonstrate that this relationship 
can be used to estimate the mass of any GBHs by simply 
scaling the Mbh value from a suitable reference source. 
We have started testing the extension of this method to 
SMBHs and the encouraging results will be published 
elsewhere (Gliozzi et al. 2010 in preparation). Here we 
apply this method to PKS 0558-504 after describing the 
basic characteristics of the BMC model and of this scal- 
ing technique. 

Although it was historically developed to describe the 
Comptonization of t hermal seed photons by a relativis- 
tic converging fiow (jTitarchuk et al.lllQQTt) . the BMC 
model is a generic Comptonization model able to de- 
scribe equally well the thermal Comptonization (i.e., the 
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inverse Compton scattering produced by electrons with 
a Maxwellian energy distribution) and the bulk motion 
Comptonization (where the seed photons are scattered 
ofi' electrons with bulk relativistic motion). The BMC 
model is characterized by 4 free parameters: 1) the tem- 
perature of the thermal seed photons fcT, 2) the energy 
spectral index a (which is related to the photon index by 
the relation F = 1 -I- a), 3) a parameter log(^) related to 
the Comptonization fraction / (.i.e., the ratio between 
the number of Compton scattered photons and the num- 
ber of seed photons) by the relation / = A/(l -t- ^), and 
4) the normalization A'bmc- 

In simple terms, the necessary steps to derive ilfsH 
with this method can be summarized as follows: 

1) Construct a F — TVbmc plot for a GBH of known mass 
and distance, which will be used as reference (hereafter 
denoted by the subscript r). 

2) Compute the normalization ratio between the target 
of interest (hereafter denoted by the subscript t) and the 
reference object A^BMCt/A^BMCr at the value of F mea- 
sured for the AGN. 

3) Derive the unknown black hole mass using the follow- 
ing equation 

MBH,t = MBH.r X (iVsMC.t/iVBMCr) X {dtjdrf X /g (1) 

where Msn.r is the black hole mass of the GBH reference 
object, A^BMC.t and A^BMCr are the respective BMC nor- 
malizations for target and reference objects, dt and dr 
are the corresponding distances, and /g — cos Or I cos Bt 
is a geometrical factor that depends on the respective 
inclination angles. 

The above formula is readily obtained by considering 
that (a) the normalization is a function of luminosity and 
distance: iVsMC oc L/d"^, and (b) the luminosity of an 
accreting BH system can be expressed by L ex r^MBH'Tt, 
where 77 is the radiative efficiency. The only assumptions 
are that different sources in the same spectral state have 
similar values of 77 and m, and that the photon index 
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is a reliable indicator for the spectral state. This is in 
broad agreement with the positive correlation between 
pho ton index and accretion rate found in brigh t AGNs 
fe.g lShemmer et al.|[200l iPapadaki s et al.ll2009[) . 

Figure [3] shows the T — Nbmc diagram for PKS 0558- 
504 in comparison with the primary reference source 
GRO J1655-40, a well known microquasar whose pa- 
rameters are tightly constrained: M bh/Mp) = 6.3 ± 0.3 , 
i = 70° ± 1°, d = 3.2 ± 0.2 kpc fGre ene et all I2OOII : 
IHiellming fc Rupenl[l995l) . For PKS 0558-504 we fitted 
the 2-10 keV energy band of the 7 relatively long seg- 
ments (with exposures up to 20 ks) that best describe 
the spectral evolution of the sou rce during the September 
2008 XMM-Newton campaign (IPapadakis et all [20T0al) 
and the 10 short segments (2 ks exposures) simultane- 
ous to the SWIFT observation, which will be described 
below. For the spectral fits, the temperature of the ther- 
mal seed photons was kept frozen at the best fit value 
obtained from the broadband SED fit (fcT = 8 - 23 eV; 
see §5.2), whereas the other parameters were free to vary. 
The 2-10 keV spectra of all segments were adequately fit- 
ted (typically Xrod ^ 0.9 - 0.95) with one BMC model 
absorbed by Galactic A^h- 

As expected, due to the longer timescales associated 
with SMBHs, the PKS 0558-504 trend in the F - Nbmc 
diagram is restricted to a small portion of the trend 
shown by GRO J1655-40 during its 2005 spectral evo- 
lution between LS and HS. Interestingly, the short-term 
spectral behavior of PKS 0558-504 appears to be consis- 
tent with GRO J1655-40 in its highly accreting state. 

Substituting in Eq. (1) a distance of 624 Mpc and 
A^BMC,t/A^BMC,r = 8 X 10~^ obtained from Fig.[3l we de- 
rive for PKS 0558-504 Mbh ^ /g 3 x 10* Mq. Although 
the inclination angle of PKS 0558-504 is unknown, the 
nearly symmetric position of the lobe-like radio struc- 
tures suggests that the system is not seen pole-on. If we 
conservatively assume i = 30° — 45°, then fa — 0.4 — 0.5, 
leading to a mass estimate of Mbh — 1-5 x 10* Mq. 

In summary, all the optically-independent methods 
consistently yield Mbh values of the order of few units 
in 10* Mq, which is about a factor ^5 larger than the 
value derived from the virial theorem. Interestingly, this 
is in agreemen t with the corrective factor proposed by 
iMarconi et al.l (|2008[ ). However a systematic study of 
the Mbh distribution of well-defined samples of NLSls 
and BLSls obtained with optically-independent meth- 
ods is necessary before drawing any general conclusion. 
Considering all four methods, Mbh in PKS 0558-504 
ranges between 6 x 10^ — 3 x 10* Af© with a mean 
of 1.8 X 10* Mq. 

5. A MULTI- WAVELENGTH VIEW OF PKS 0558-504 

5.1. Short-term variability 

Before investigating the properties of SEDs derived 
from the multi- wavelength campaign in September 2008, 
it is important to assess the presence of variability at the 
different wavelengths probed by SWIFT. The light curves 
obtained from the six UVOT filters and from the XRT 
are shown in Fig. S) The plotted optical and UV time se- 
ries are the fiuxes in units of 10"^** erg cm ~^s~^A~^ and 
are c orrected for Galactic absorption (see iRoming et al.l 
l2009l for details). To be conservative, we have used error- 
bars of 3cr to assess the presence of significant variabil- 



ity with a x^ test (the same conclusions are obtained 
with 2a errors), because the Icr statistical errors undcr- 
estir nate the true unc ertainty on the optical-UV fluxes 
feee lPoole et al.ll2008 ^. In the top panel of Fig. g] in ad- 
dition to the XRT light curve, we show the occurrence 
and duration of the XMM-Newton EPIC pn observations. 
Sub-intervals of the duration of 2 ks simultaneous to the 
S'WTFT' pointings were used for the SED analysis. 

According to a x^ test no signiflcant variability is 
present in the U, B, V bands {Py,2 ~ 0.9 — 0.95) and 
only marginal variability is detected in the UVWl band 
(P^2 ~ 0.2), as suggested by a visual inspection of the 
four bottom panels of Fig. 21 On the other hand, the 
two higher energy UV bands (UVM2 and UVW2) and 
the X-ray range all show statistically significant vari- 
ability (P-^2 < lO"**). There are however qualitative 
and quantitative differences between the UV and X-ray 
light curves: the former show only moderate variabil- 

ity, Fvar.M2 = (3.0 ± 1.4)%, Fvar,W2 = (4.1 ± 0.5)% 

(where -FVar is the fractional variability that measures 
the normalized variance corrected for the statistical un- 
certainties), which appears to be associated with only 
two data points (specifically the third and fourth points 
in Fig.S]) that are significantly lower than the others. On 
the other hand, the XRT light curve shows strong vari- 
ability (Fvar,XRT = (28.6 ±0.9)%) throughout the entire 
8-day interval. 

In conclusion, the 5* WffT campaign in September 2008 
confirms the presence of strong variability in the X-ray 
band of PKS 0558-504 and reveals the presence of mod- 
erate variability in the UV bands, whereas the optical 
bands appear to be consistent with the hypothesis of con- 
stant flux over short timescales. 

5.2. Broad-band SED 

One of the best ways to get insights into the energetic 
of the central engine in PKS 0558-504 and in AGNs in 
general is to study the broadband SED, as de monstrated 
by nu merous successful past studies (e.g. see lElvis et al.l 
|1994|) . However, one of the major problems in this fleld 
has been the lack of truly contemporaneous observations 
in different energy bands, which could lead to inaccurate 
results in virtue of the fast variability of AGNs. Re- 
cently this kind of studies has markedly improved with 
the advent of SWIFT thanks to its ability of observing 
simultaneously in several energy bands coupled with its 
highly flexible schedule. 

Here we analyze the contemporaneous SEDs of 
PKS 0558-504 obtained with the SWIFT UVOT and 
XMM-Newton EPIC pn during a simultaneous campaign 
in September 2008. The choice of using 2000 s segments 
of the EPIC pn instead of the SWIFT XRT data is dic- 
tated by the much higher throughput of the EPIC camera 
that allows a more detailed spectral analysis of the X-ray 
spectrum. 

Following the pro cedure adopted by 

IVasudcvan fc FabianI (|2009[ ). we converted the UVOT 
fluxes using the FLX2XSP command in FTOOLS. The 
resulting spectra were then combined with the 0.3-10 
keV EPIC spectra and fitted in XSPEC with a model that 
comprises a disk and two Comptonization components: 
DISKPN+WABS(BMC-hBMC). Only the Comptoniza- 
tion models are absorbed by the Galactic A^h, because 
the optical and UV data were already corrected for 
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absorption. 

As explained in §4 the BMC is a simple but compre- 
hensive Comptonization model that can fit both thermal 
and bulk Comptonization processes. The use of a Comp- 
tonization model to parametrize the soft excess below 2 
keV is guided by a thorough spectral analysis of the high- 
est quality XMM-Newton spectra, which suggests that 
this component may arise from a hot " skin" in the inner- 
most part of the disk in objects that accrete at rates clos e 
or above the Eddington limit (jPapadakis et al.ll2010aD . 
For the hard X-rays we also used a BMC model instead of 
the phenomenological power law model (PL) because the 
BMC parameters are computed in a self-consistent way, 
and, unlike the PL, the power law produced by BMC 
does not extend to arbitrarily low energies and thus does 
not affect the normalization of the thermal component 
nor the amount of local absorpt ion. 

The DISKPN model (Gierlihski et alJ[l999t ) has three 
parameters: Tmax, ^in, and the normalization that de- 
pends on the black hole mass, the distance, the inclina- 



tion angle i, and the color factor /3. The latter two quan- 
tities, which are unknown for PKS 0558-504, were kept 
fixed to i = 0° and /3 = 1, since they only marginally af- 
fect the luminosity derived from the direct integration of 
the SED (see Vasud evan & F abian 2009) . In order to de- 
termine the maximum temperature of the accretion disk, 
we first fitted the UVOT data only with the DISKPN 
model, fixing i?in = 6 Rq and the normalization to two 
different values corresponding to Mbh = 6 x 10^ Mq and 
Mbh — 2.5 X 10^ Mq, respectively. The resulting best 
fit temperatures were respectively ^ 23 eV and ~ 8 eV. 
If the normalization of the DISKPN model is left free 
to vary, the resulting value derived for the BH mass is 
A/bh — 2.5 — 2. 7x 10* Mq, which in good agreement with 
the values inferred from the three model-independent 
techniques described in Sect. 4. 

We then fitted the X-ray spectra by linking the tem- 
perature of seed photons of both Comptonization models 
to Tmax- The resulting best fit parameters for the BMC 
components and the corresponding luminosities are re- 
ported in Table 3 for the conservative assumption that 
Mbh = 2.5 x 10* Mq. The values of the reduced x^ for 
the X-ray spectral fits range between 0.9 and 1.2, whereas 
for the overall broadband fits Xr are of the order of 1.4- 
1.5. Note that the BMC parameters and the X-ray lumi- 
nosity are almost insensitive of the choice of fcT, whereas 
the bolometric luminosity Lboi — ^o.oi-ioo keV is signifi- 
cantly affected by this choice: Lboi increases by a factor 
of -3 when kT = 23 eV (i.e., for Mrh = 6 x 10^ Mq). 

Since a very detailed X-ray spectral analysis has 
been already performe d and reported elsewhere (see 
iPapadakis et al.l l2010a[ ). here we focus on the broad- 
band SED. Figure [5] illustrates three different SEDs cor- 
responding respectively to (a) segment 001 in the top 
panel, which represents the case of average optical-UV 
fiux and average X-ray emission; (b) segment 002 in the 
middle panel, which refers the highest X-ray emission 
level; and (c) segment 003 in the bottom panel that rep- 
resents the case of low optical/UV emission. Whereas 
the X-ray emission varies considerably from day to day 
with luminosity changes up to a factor of ^2, the changes 
in the optical/UV band are barely noticeable and corre- 
spond to variations of Lboi lower than 10%. 

In conclusion, the overall SEDs of PKS 0558-504 dur- 
ing the September 2008 campaign are reasonably well 
fitted with a disk model parameterizing the optical/UV 
emission and two Comptonization components to de- 
scribe the X-ray spectra. The SED appears to be domi- 
nated by a nearly constant disk component, with a vari- 
able X-ray contribution which represents about 1/100 of 
the Lboi with the assumption that Mbh = 2.5 x 10* Mq. 
The X-ray contribution is further reduced by a factor of 
-3if A/bh =6 X IO^Mq. 

6. DISCUSSION 

6.1. Jet contribution 

The main goal of this work is to investigate the en- 
ergetic of the radio-loud NLSl PKS 0558-504 and shed 
some light on the nature of its central engine. To this 
end, it is first necessary to assess the role played by the 
jet and more specifically whether the jet may dominate 
the high energy emission of PKS 0558-504. The findings 
based on the ATCA observations reported in §3 - the 
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TABLE 3 

X-RAY SPECTRAL PARAMETERS AND LUMINOSITY OF PKS 0558-504 
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logAs 
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Columns Table 3: 1— Segment of the SWIFT campaign. 2— Spectral index of the BMC model fitting the soft excess belcw 2 IceV. 3— Measure 
of the Comptonization fraction (/ — a/(l + A)) for the soft excess. 4— Spectral index of the BMC model for the 2-10 keV range. 5— Measure of 
the Comptonization fraction for the 2-10 keV range. 6— Soft X-ray luminosity. 7— Hard X-ray luminosity. 8— Boloinetric luminosity 
Notes: The bolometric luminosity is obtained by assuming A/bh — 2.5 X 10 Mq, which leads to a best fit temperature of kT 2± 8 eV. If 
Mbh — 6 X 10 Mq, then kT ^ 23 eV and the bolometric luminosity is a factor ~ 3 higher than the values reported in the Table. 



radio-loudness measurements consistent with radio-quiet 
Seyfert galaxies and the extended and symmetric FRI- 
hke radio structure which is at odds with the compact 
radio emission typical of very radio- loud NLSls- disfa- 
vor a jet-dominated scenario, which instead appears to 
be the most likely scenario for the very radio-loud NLSls 
recently detected with the Fermi/LAT. 

This preliminary conclusion is supported by several 
independent pieces of evidence from studies at higher 
energies, which can be summarized as follows: (1) the 
X-ray spectral and temporal variability appears to be 
fully consistent with th e typical behavior of radio-quiet 
Seyfert galaxies (jPapad akis ct al. 2010a,b); (2) the long- 
term spectral variability in the 2-15 keV energy band 
is inconsistent with th e characterist i c tre nd observed 
in blazars with RXTE (jGliozzi et all 120071 ): (3) unlike 
highly beamed AGNs, PKS 0558-504 is not detected at 
very high energy neither in the GeV range (the first year 
Fermi LAT catalog) nor at TeV energies (HESS; Giebels 
private communication). Combining all the information 
available for PKS 0558-504 it is thus reasonable to con- 
clude that the jet emission does not play a relevant role 
beyond the radio band. Therefore, for the rest of the pa- 
per we will assume that the bulk of the optical-UV- X-ray 
emission is associate with the accretion flow only. 

6.2. Energetic of PKS 0558-504 

The main result from the analysis of the broad- 
band SED (described in §5) is that the emission of 
PKS 0558-504 is dominated by a weakly variable UV 
bump (i^var,uv — 4%), which is commonly associated 
with the emission from the accretion disk. Conversely, 
the X-ray radiation, which we have interpreted as emis- 
sion from two Comptonization components, appears to 
be highly variable (Fvar.x — 30%) but encompasses only 
a small fraction of the total emission. It is instructive to 
assess in a more quantitative way which of the variable 
components (UV or X-rays) is energetically more impor- 



tant. Since -FVar provides a measurement of the variabil- 
ity normalized over the average flux, a measurement of 
the "true" variability (hereafter varuy and varx, respec- 
tively) can be obtained by multiplying i^Var by the respec- 
tive average flux. This test yields varx/varuy = 1.3, in- 
dicating that the X-ray component is the dominant vari- 
able component. 

In order to get helpful insights into the nature of the 
central engine in PKS 0558-504 it is crucial to constrain 
its accretion rate in term of Eddington units. This in- 
formation can be readily obtained by knowing the mass 
of the central black hole and measuring the bolomet- 
ric luminosity. Depending on the method used (see §4) 
the black hole mass of PKS 0558-504 may range be- 
tween 6 X 10^ M0 and 3 x 10^ M©. To be conserva- 
tive, in the rest of the paper we will assume A/bh = 
2.5 X 10* Mq which corresponds to an Eddington lumi- 
nosity of ^Edd — 3.25 X 10''^ erg s^^ (note that this 
choice of Mbh does not affect the main findings nor 
our conclusions). Combining this value with the aver- 
age bolometric luminosity derived from the direct inte- 
gration of the broadband SED (see Table 3) leads to an 
accretion rate of Asdd = ^boi/^Edd — 1-7 that strongly 
suggests that PKS 0558-504 accretes at super-Eddington 
rate. Note that with the less conservative assumption of 
Mbh = 6 X 10^ M© (which corresponds to a lower Ed- 
dington luminosity and a higher bolometric luminosity), 
the accretion rate in Eddington units would be of the 
order of AEdd — 25. 

Having measured the X-ray and bolometric luminos- 
ity over a period of ten days, it is possible to compute 
not only the average bolometric correction K2-10 keV = 
^2-10 kev/iboi but also to investigate if it varies over 
time. From Table 3 we infer that K2-10 keV indeed varies 
by a factor of '•^2.5 ranging between 108 and 249, with 
an average value of 162 ± 39 (where the error quoted 
is la). The large value of K2-10 keV appears to be in 
qualitative agreement with recent studies that suggest 
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Fig. 5. — Deconvolved EPIC pn spectra in the 0.4-10 keV energy 
range, combined with UVOT data. The overall SEDs are fitted with 
a diskpn model to parametrize the optical-UV emission and two BMC 
models to characterize the X-ray spectrum. The BMC models are ab- 
sorbed by Galactic A^hi whereas the optical/UV data have been already 
corrected for Galactic absorption. The spectrum in the top panel (a) 
represents the case of average optical-UV flux and average X-ray emis- 
sion. The middle panel (b) illustrates the case of high X-ray emission, 
whereas the bottom panel (c) represents the case of relatively low op- 
tical/UV emission. 



that the X-ray bolometric correcti on is directly propor- 
tiona l to the Eddington ratio (e.g. IVasudevan fc FabianI 
120071 12009). A more quantitative compariso n can be per- 
formed using the recent results presented by iLusso et al.l 
(|2010) . which are based on the optical and X-ray study 
of a sample of 545 X-ray selected type 1 AGNs from 
the XMM-COSMOS survey. More specifically, plugging 
^Edd — 1-7 (the value of the Eddington r atio derived 
for PKS 0558-504) into equation (14) from iLusso et al.l 
(|2010i) we derive k.2-w keV = 151^20' which is fully con- 
sistent with the average value empirically derived from 
the SED. 

Since PKS 0558-504 possesses a radio jet, it is instruc- 
tive to estimate the the inflow accretion rate and try to 
compare it to the the outflow mass rate associated with 
the jet, and similarly compare the accretion luminosity 
to the kinetic power. With the reasonable assumption 
that Lboi = iacer = TlMi^c^ , we derive an accretion in- 
flow rate of Afjn c± 1 M0/yr, assuming a typical radiative 
efficiency of 77 = 0.1 (larger values of M are obtained if 
photon trapping effects play a signifi cant role causing 

2 < to d ecrease considerably). Following iMerloni fc Heina 
2007( ) it is possible to derive an estimate of the kinetic 
power of the jet for radio galaxies by making use of the 
formula logLkin — 0.81 log Lr -|- 11.9, where Lr is the 
nuclear radio luminosity at 5 GHz. Using the value mea- 
sured with the ATCA, Lsghz = 1-7 x 10^^ erg s"\ we 
derive a kinetic power of the order of Lkin — 2 x 10^^ erg 
s~^. The latter is considerably lower than the luminosity 
associated with the accretion flow, which is of the order 
of 5 X lO**^ erg s~^. We can estimate the outflow mass 
rate from Lkin = (r— l)MoutC^ by assuming an appropri- 
ate value for the bulk Lorentz factor. With F = 3, which 
is a typical value for radio galaxies and consistent with 
the value inferred from the brightness temperature de- 
rived from the radio variability (Kedziora-Chudczer et al. 
2010, in preparation), we derive Mout — 2 x 10~^ MQ/yr, 
which is about two orders of magnitude lower than the 
value inferred for the accretion rate Mjn. 

6.3. Analogy with GBHs 

Our recent study based on a 1-year RXTE monitoring 
campaign suggests that the spectral variability behavior 
of PKS 0558-504 in the 2-15 keV energy band mimics 
GBHs in their highly accreting IS f Gliozzi et all 12(307). 
We can now test this analogy in a more quantitative way. 
First, we can compare the results from the energetic sec- 
ti on relative to the dis k-jet interaction with the findings 
of lOhsuga et al.l (|2009l ) that are derived from radiation- 
magneto-hydrodynamic simulations of a stellar BH sys- 
tem accreting at the Eddington rate. The general agree- 
ment with the simulations (specifically, Afout '^ -^^^in and 
ikin < -^accr) Icnds further support to the hypothesis 
that PKS 0558-504 is a large scale analog of a GBH in 
a highly accreting state. 

The analogy with the IS/VHS is supported by the di- 
rect comparison of PKS 0558-504 with GRO J1655-40 
in the F — A^bmc diagram that we used to estimate the 
black hole mass (see Fig. [3]). Specifically, the common 
spectral variability trend occurs when GRO J1655-40 is 
switching from the LS to the HS (i.e., the source is indeed 
in the IS) duri ng the rise of the outburst in March 2005 
(see details in IShaposhnikov et al] I2007D . Interestingly, 
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the X-ray outburst of GRO J1655-40 was accompanied 
by enhanced emission in the radio band that faded before 
the peak of the X-ray outburs t . Bas ed on the values re- 
ported by IShaposhnikov et al.l ()20 07^. the ratio between 
radio/X-ray and radio/optical fluxes places GRO J1655- 
40 in the same lower-left quadrant of the radio-loudness 
plot shown in Fig. [U for PKS 0558-504. 

Further support to this conclusion comes from the 
comparison of the broadband spectral index aox of 
PKS 0558-504 with the correspondin g quantity anRHi 
which was recently introduced by ISobolewska et al.l 
()2009| ) to characterize the broadband spectrum of GBHs 
and allows a direct comparison with AGNs. During our 
short multi-wavelength campaign in October 2008, aox 
ranged between 1.25-1.4 while the photon index varied 
from 2.1-2.25. Importantly, this is the very range corre- 
spo nding to the interm ediate state in GBHs (see Figure 
6 in ISobolewska et al.l [20091 ). 

Finally, it is instructive to compare PKS 0558-504 to 
GRS 1915-fl05, which is thought to be the only mi- 
croquasar that regularly accretes at or above the Ed- 
dington limit. Interestingly, during its spectral transi- 
tion hard-to-soft the X-ray spectra of GRS 1915-1-105 are 
well fitted by two BMC components (a hard component 
with Fi =1.7-3 and a soft component with Fi =2.7- 
4.2) which evolve followin g the trend shown in Fig. |3| 
(jTitarchuk fc Seifinall2009l ). Also for GRS 1915-^105 the 
location of the IS/VHS in the F — Nbmc diagram is con- 
sistent with the position of PKS 0558-504 (see Figure 13 
of .Titarchuk fc Sci fina 2009). 

6.4. Conclusion 

In conclusion, our multiwavelength analysis of the 
emission properties of PKS 0558-504 ranging between 
the radio and the hard X-rays reveals that this source 
has a radio jet, yet it is accreting at the Eddington level 
or above. All our findings appear to independently lend 
support to the hypothesis that PKS 0558-504 may be a 
large scale analog of the IS observed in GBHs. But can 
we extend this analogy to the whole class of NLSls? In 
other words, can PKS 0558-504 be considered a proto- 
type of NLSls? At first sight, the "radio-loudness" of 
PKS 0558-504 seems to set this source in a small sub- 
class of radio-loud NLSl. However, _R = 10 does not 
represent a strict boundary for the radio-loudness; in- 
deed the radio-loudness of PKS 0558-504 appears to be 
fully consistent with that of radio-quiet Seyfert galaxies 
according to both R and Rx (see Fig. [2]). Bearing in mind 
that PKS 0558-504 is by far the brightest NLSl in the 
X-ray band and that its X-ray luminosity is nearly 100 
times the typical value of "normal" radio-quiet NLSls, 



similar values of the order of i?^; < lO"'^ for "normal" 
NLSls would correspond to radio fluxes of the order of 1 
mjy or less. This would explain why only a small minor- 
ity of NLSls is detected in the radio band. Additionally, 
it must be pointed out that in the IS the radio ejections 
are transients, therefore, genuine radio-quiet NLSl might 
still represent large scale analogs of the IS. 

It is also interesting to compare PKS 0558-504 with the 
radio-loud NLSls recently detected at 7-rays. The latter 
are characterized by large values of the radio-loudness 
{R > 50), by compact radio emission, and by SEDs 
that are dominated by two broad bumps typical of jet- 
dominated sources. Indeed, it has been argued that these 
very radio-loud NLSls galaxies represent a the third class 
of 7-ray emitting AGNs along with blazars and radio 
galaxies. In this perspective, PKS 0558-504 can play 
a crucial role since it may represent the prototype of 
NLSls with a jet seen at large viewing angles. In other 
words, PKS 0558-504 may be the the first representa- 
tive of non-jet-dominated radio-loud NLSls which are the 
parent population of the very radio-loud NLSls detected 
by the Fermi/LAT, similarly to radio galaxies that are 
considered the non-beamed parent population of blazars. 

In summary, we can hypothesize that PKS 0558-504 
is indeed a prototypical NLSl. Since the jet appears to 
dominate only in the radio band, the emission associ- 
ated with the accretion fiow can be cleanly separated by 
that related to the jet. This makes PKS 0558-504 an 
ideal system to study the disk-jet interaction. In this 
framework, important insights into the physical condi- 
tions leading to the formation of jets in highly accreting 
BH systems can be obtained by a systematic comparison 
of the spectral and temporal properties of PKS 0558-504 
with the corresponding properties of other highly accret- 
ing NLSl that do not show any evidence of radio jets as 
well as with NLSls that are jet-dominated. We plan to 
pursue this approach in our future work. 
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